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Abstract. It is well known that the Husimi Q-function of the signal field can actually be 
measured by the eight-port homodyne detection technique, provided that the reference beam 
(used for homodyne detection) is a very strong coherent field so that it can be treated classically, 
see e.g. [Hj. Using recent rigorous results on the quantum theory of homodyne detection 
observables |12], we show that any phase space observable, and not only the Q-function, can be 
obtained as a high amplitude limit of the signal observable actually measured by an eight-port 
homodyne detector. The proof of this fact does not involve any classicality assumption. 



1. Introduction 

Covariant phase space observables, as positive operator measures, play an important role in 
the foundations of quantum mechanics. In particular, their importance for the approximate 
joint measurements of position and momentum observables has long been recognizee^ and it 
has recently been shown [T7| that for any approximate joint measurement of position and mo- 
mentum of a quantum object there is a covariant phase space observable with improved degrees 
of approximations Jl The mathematical structure of the covariant phase space observables is 
also completely known: they correspond one-to-one onto the positive operators of trace one 
(acting on the Hilbert space of the quantum object in question), and they have an operator 
density defined by the Weyl operators and the positive trace-one operator in question, see eq. 
@ belowj^ Moreover, the realization of the covariant phase space observables (associated with 
one-dimensional projection operators) as sequential position-momentum measurements, with 
the first measurement as an approximately repeatable measurement, has recently been demon- 
strated in [7] following the pioneering work of Davies What remains then is the question 
of the experimental implementation of the phase space observables. 

It is well-known that using the classical approximation of high amplitude reference beam, the 
Husimi Q-function (i.e. the phase space observable generated by the vacuum state operator) is 
obtained as the measured observable for the so called eight-port homodyne detector (see fTS] , 
and [m p. 147-155]). In this paper we show, using the recent results [I2] on the balanced 
homodyne detection observables, that any covariant phase space observable can be obtained, in 
a mathematically rigorous sense, as the high amplitude limit of signal observables determined 
by the eight-port homodyne detection scheme. 
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Figure 1. The eight-port homodyne detector 

2. The eight-port homodyne detector 

The detector involves four modes as indicated in the picture, and we will denote the associated 
(complex separable) Hilbert spaces accordingly by Hi, H2, H3, Tt^. Mode 1 corresponds to 
the signal field (i.e. the object system with respect to which the measured observable will be 
interpreted), the input state for mode 2 serves as a parameter which determines (as will be 
seen below) the phase space observable to be measured, and mode 4 is the reference beam in a 
coherent state. (The input for mode 3 is left empty, corresponding to the vacuum state.) 

We fix the photon number bases {|n) | G N} for each Hi, so that the annihilation oper- 
ators ttj, as well as the quadratures Qj = -^{ci*j + aj), Pj = — aj) and photon number 
operators Nj = a*aj are defined for each mode j = 1,2,3,4. The bar above denotes the clo- 
sure of an operator. We will also sometimes use the coordinate representation for each Hj (i.e 
Hj ~ L'^iM.), with In) associated with nth Hermite function). Then Qj and Pj act as usual 
position and momentum operators: {Qj'ilj){x) = x^p^x), and {Pjifj){x) = i/) G L^(M). 
For any selfadjoint operator A on a Hilbert space, we let P^ denote its spectral measure. 



^See, for instance, the monographs [8l flOl [2t [16] . 

^For details of these concepts as well as for a further analysis of these results, see the above quoted work of 
Werner [17] as well as the subsequent developments [Sl [HI . 

^This result is due to Holevo [TI] and Werner [18], alternative proofs with different techniques were recently 
given in [5] and |13) . 
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The photon detectors Dj shown in the picture are considered to be ideal, so that each detector 
Dj measures the sharp photon number Nj. The phase shifter in mode 4 is represented by the 
unitary operator e*'^^'', where (j) is the shift. 

There are four 50-50-beam splitters S12, -B43, f/13, f/24, each of which is defined by its acting 
in the coordinate representation (see e.g. [T5]): 

L\m?) 9 ^ ^ ( (xi, X2) ^(71(^1 + ^2), ^(-2:1 + X2)) ) G L2(M2), 

Under this transform, the coherent states change according to 

(1) l«) ® I/?) ^I7l(«-/9))®l7l(« + /?))• 

In the picture, the dashed line in each beam splitter indicates the input port of the "primary 
mode", i.e. the mode associated with the first component of the tensor product L^(]R) (S)L^(M) ~ 
L^(M^) in the above description. The beam splitters S12, ^43, f/13 and t/24 are indexed so that 
the first index indicates the primary mode. (For instance, i?43 acts such that an input two-mode 
state ® e 7-^3 ® 7-^4 is interpreted as \l/(zi,X2) := il){xi)^p{x2) in the above definition.) 

The unitary operator describing the entire transform caused by the combination of the beam 
splitters together with the phase shifter is 

(^13®t^24)e*'^'^^^®^^(5i2®S43). 

Here it is understood that 

7^1 ® 7^2 ® ^3 ® ^4 ^ (^1 ® ^2) ® (^3 ® ^4) ^ {'Hi ® Hs) ® (Hi ® Hi), 

and we will freely use these isometrics, without explicit indication, when the ordering is clear 
from the context. 

Let I a/22;) be the coherent input state for mode 4. We choose to detect the scaled number 
differences |f|^i3 and j7|-?V24i where 

:= h^Ns- Ni0 h, N^^ := h^N^- h, 
so that the detection statistics are described by the biobservable 

{X,Y) ^ pl^l"'^r3(X) (g) p\^\~'^2,(Y) 

acting on the entire four-mode field. 

Let the input states of modes 1 and 2 be T and S, respectively, so that the input for the 
four-mode field is T ® ® |0)(0| ® \^z) i^\plz\ (with the natural ordering of tensor products). 
The action of 543, as well as the phase shifter can be calculated explicitly in terms of coherent 
states: 

(/3 ® e'^^^)S43(|0) ® \^Piz)) = \z) ® \e''t'z). 

(Notice that here mode 4 is the primary mode, so that equation ^ is applied with the tensor 
product order reversed.) Now the state of the system after the transform can be written as 

where 

Wt,s,z,<i> := 5i2(T® 5)5*2 ® \z){z\ \ze''^){ze% 
(Here the tensor product is written in the original order.) 
The detection statistics are given by the probability bimeasures 

p'/'%X,Y) := Tr[t/i3 ® U2aWt,s,z,^UI^ ® t/*4Pl^l"'^»(X) ® pl^l^'^-(F)]. 
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Since -^U^^N^sUrs = ^13, and -^UliN^JJ2A = A24, where e.g. A13 = ^(ai a| + ® ag), 
we have simply 

Let : i3(M^) LiHi ®l-L2®'Hz® 'Hi) denote the unique spectral measure extending the 
biobservable (X, F) ^ Pl"l"'^"(X) ® pl^l"'^24(y)^ go ^hat 

p'/^f{X,Y) = Tr[W^^,5,.,0P^(^(X X Y))] 

for any X,Y e B{R). 

In our measurement, the input state S of mode 2, as well as the coherent input state \\/2z) 
of mode 4, and the phase shift (f), are regarded as fixed parameters, while T is the state of the 
object system, i.e. mode 1. Accordingly, we define an observable G^'^'"^ : i3(M^) L{T-Ci) via 

(2) Tr[rG^'^'^Z)] = Ti[WT,s,z,<pP'ij^Z)l Z G ^(M^). 

(Since TT[WT,s,z,^P%j^Z)] = Tt[T 0S® \z){z\ ® \ze^'^){ze'^\Bl2 ® 134/^^(75^)512 ® 134], it 
follows by standard duality and convergence arguments that the observable G^'^''^ exists and is 
uniquely determined by the above formula.) The observable measured by the eight port 
homodyne detector is thus G^ '^'^. 

In order to consider the limit \z\ 00, we need to express G^'^'"^ in terms of the single 
homodyne detector observables and Ef\ where Ef{X) := {V}^)*P\'\~'^^^{X)V^^^, V}^ : 
Hi ^ Hi ® H3 is the isometry 1— > ® jz), and E^'^' (Y) is defined in an analogous way by 
using yl24 and Kfgf^ (see [l2]). We let S{Hi) denote the set of states (that is, positive operators 
of trace one) of the mode described by the Hilbert space Hi, i = 1,2, 3, 4. 

Lemma 1. For any T G S{Hi), X,Y & 'B(ffi), and an S E S{H2), we have 

(3) Tr[rG^'^'^(X X Y)] = Tr[Si2(r ® S)Bl2El{^^X) ® Ef'{j^Y)]. 

Proof. Let X,Y E B{R). Assume first that T = P[(p] and S = P[ip] for some unit vectors 
G Hi, V e H2. Write 

00 

n,m=0 

where {r^^ | ?7, G N} is an orthonormal basis for i = 1, 2. This series converges in the norm 
of Hi ® 7^2, so by using the definition we get 

Tr[TG'"'^''^(V2(A: x Y))] = (Pi2(^ ® ^A) ® \z) ® \e''^z) P'{X x Y)Bi2{^ ® ^) ® 1^) ® \e''^z) 
00 

^^n'm'^ri ®rt® \z) ® |ze*'^)|pl^l"^^='(X) ® p\^\'"^^\Y)ji, ® ril, ® \z) ® W^) 

n,m,n' ,m'=0 
00 

= E ^c^'rn'{ril\{V^yp\'\''^^^{X)V^'vl) 



n,m,n' ,m'=0 
00 



n,m,n' ,m'=0 



{Bui^ ® ^)l^i(^) ® Ef\Y)Bi2{v ® ^)) = Tr[Pi2(T ® 5)Pi*2^r(X) ® ^2''"'(r)]. 
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As for the general case, note that both sides of (JS]) are trace norm continuous when regarded as 
functions of T (for fixed 5") , and the same is true when they are regarded as functions of S for 
fixed T. Hence, by using the spectral resolutions for S and T (which are trace norm convergent 
weighted sums of one-dimensional projectors), one establishes the claim. □ 

3. The high amplitude limit 

Now we proceed to describe the limit \z\ oo, in the case where z = r > 0, and = f • In 
order to simplify the notation, we will drop the subscript for mode 1 operators (e.g. Q := Qi) 
from now on. We need the following general definition, which we introduced in |12| . 



Definition 1. Let H he a Hilbert space, fl a metric space with Borel cr-algebra B{Q), and 
E"- : B{Q) LiTi) a semispectral measure for each n. G N. We say that the sequence (i?") 
converges to a semispectral measure E : B{Vt) — > L(l-L) weakly in the sense of probabilities, if 

lim E'^iX) = E{X) 

n— ►oo 

in the weak operator topology, for all X G B{VL) such that E{dX) = 0, where dX is the 
boundary of the set X. (Recall that dX is the intersection of the closures of X and its 
complement.) 

We used this definition in the context where = M, and proved that for any 6 G [0,27r), 
and any sequence (r„) of positive numbers converging to infinity, the sequence (i?''"'^"^)„gN of 
homodyne detector observables converges to P'^" weakly in the sense of probabilities, where 
Qg = e^^^Qe~^^^ is the rotated quadrature [l2l P- 17]. Since the spectral measure of each 
quadrature is absolutely continuous with the Lebesgue measure A, the condition P^^^dX) = 
for an X G B{R) is reduced to the condition \{dX) = 0. Now we need only the particular 
choices = and = f , when the limiting observables are the spectral measures of Q and P, 
respectively. 

For any positive operator S G LiTii) of trace one, let E^ : i3(R^) L{Hi) denote the phase 
space observable generated by 5, i.e. 

(4) ES{Z) = ^ [ W{q,p)SW{q,prdqdp, 

^TT Jz 

where W{q,p) = e*5'?Pe~*^^e*P'3, q,p e R are the Weyl operators. Let C : H2 ^ Hi denote 
the (antiunitary) conjugation operator i/j (x \—>- i/^{x))) Notice that here C is interpreted a 
map from H2 — L^(M) to Hi ~ L^(R). Since C is antiunitary, it follows that for any positive 
operator 5* G L{H2) of trace one, the map CSC~^ is a bounded linear positive operator in Hi 
with unit trace. 

The following well-known result can be found e.g from [2], p. 195-196]. We reproduce it here 
so as to make sure that the notations fit together correctly. 

Lemma 2. Let S G L{H2) be a positive operator of trace one, and T G S(Hi). Then 

(5) Ti[TE^^^'\X xY)] = Ti[Bi2{T®S)Bl^P^{j^X)®P'''{j^Y)], X,Y e BiR). 

Proof. We have P2 = F~^Q2F, where F is the Fourier-Plancherel operator. If G L^(M) fl 
L^(R), the operator F acts as the Fourier transform: 
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The relation 



V27T 

holds for all ^ G Hi, t] G Ti.2, = 1, all p G M and almost all g G M, with the function 
X T]{x - q)^{x) belonging to L^(R) n L'^{R) for almost q eR (see e.g. [H, p. 47, 49]). 

First, let S = P[ip] for some unit vector ^ G 7^2, and let ip E Hi. Because of the above 
relation, we have 

V2^{WiV2q,V2p)Cij\ip) = V2e''iP[F[ilj{- - V2q)ip]iV2p) 



V2^= [ e-'P^^''-''^ipix)ij{x - V2g) dx 
v27r J 

e-^Py^ij=^iq + y))^ij=^{-q + y))dy 



V27r 

= F[[5i2(¥^®V^)](g,-)](p) 
for all p G M and almost all g G M, and, subsequently 

= {{I ® F)Bi2{^ ® ^)|pQ(_i.X) ® p02(-l.y)(/ ® F)5i2(v:^ ® ^)) 

|[F[Si2('^®^)](g,-)](p)l'rfgrfp 



1 

V2 



(XxY) 



2\{W{V2q,V2p)Cip\ip)\^ dqdp 



= ^ [ \{W{q,p)Cij\^)\'dqdp 

Hence, we have shown that (l5l) holds when S = Plip] and T = Since Tt[TE^'{Z)] = 

Tt[S' E'^ {— Z)] for any Z G i3(R^) and a positive operator S' of trace one, the linear maps 
T I— > Tt[TE'"^'" (Z)] and 5 i— > Tr[Ti?'-^'^'" (Z)] are trace-norm continuous. Clearly also the 
right hand side of (l5l) depends linearly and continuously on both S and T, so that the proof is 
completed by applying the spectral representations for S and T. □ 

The convergence in the high amplitude limit is characterized by the following Proposition, 
where Definition [T] is used in the case where f2 = M^. 

Proposition 1. Let S G L(H2) be any positive operator of trace one, and let (r„) be any 

r S — 

sequence of positive numbers converging to infinity. Then the sequence (G^"' "^)nm converges 
to the phase space observable E^^^ weakly in the sense of probabilities. 



Proof. Fix (r„) to be a sequence of positive numbers converging to infinity. Let X, y G 
be such that X{dX) = X{dY) = 0, where A is the Lebesgue measure of M. According to 
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the discussion at the beginning of the section, (£'["( ^X))„gN converges to P'^{^X) in the 
weak operator topology of L{Hi), and (£'2"*(^F))„eN converges to P^^{-^Y) in the weak 
operator topology of L('H2). Since the norms of all these operators are bounded by 1, it 
follows that the tensor product operator sequence i?2"*(:^l^))neN converges to 

P'^{-^X) P^'^{^Y) in the weak operator topology of Hi ^0.2- The boundedness of the 
operator norms further implies that the latter sequence actually converges ultraweakly. Hence, 
it follows from Lemma [H and [2] that for any T G 5 [Hi ) , we have 

lim Tr[rG'"-^'i(X x Y)] = Ti[Bi2{T ® S)Bl^P^{^X) ® P^%^Y)] = Tt[TE^^^~\x x Y)]. 

n—>oo 

Since the family 

{X xY \ X,Y e 13{R), X{dX) = X{dY) = 0} 

is closed under finite intersections (note that d{Xi fl X2) C d{Xi) U d{X2), Xi,X2 G i3(]R)), 
and includes a neighborhood base (for the usual topology of M^) of any point x G M^, it follows 

from [H Corollary 1, p. 14] that for any T G SiHi), the sequence {Z ^ Tr[rG''"'^'t (Z)])„gN of 
probability measures on converges weakly to the probability measure Z \—> Tr[TE^^^ (Z)]. 
(Recall that a sequence (/i„) of probability measures on is said to converge weakly to a 
probability measure fi, if lim„ J f dfin = J f dfj, for each bounded continuous real function / [H 
p. 11].) An application of Proposition 10 of [T2j now completes the proof. □ 

Remark 1. Since any phase space observable is absolutely continuous with respect to the 
Lebesgue measure, it follows that 

w - lim G^-'^'f (Z) = E^^^"\Z) 
for any Z G i3(]R^) with \^{dZ) = 0, where is the two-dimensional Lebesgue measure. 
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